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Abstract 

The growth of balhstic aggregates on deterministic fractal substrates is studied by means of 
numerical simulations. First, we attempt the description of the evolving interface of the aggregates 
by applying the well-established Family- Vicsek dynamic scaling approach. Systematic deviations 
from that standard scaling law are observed, suggesting that significant scaling corrections have to 
be introduced in order to achieve a more accurate understanding of the behavior of the interface. 
Subsequently, we study the internal structure of the growing aggregates that can be rationalized 
in terms of the scaling behavior of frozen trees, i.e., structures inhibited for further growth, lying 
below the growing interface. It is shown that the rms height {hs) and width {wg) of the trees 
of size s obey power laws of the form hs oc s'^li and Wg oc s'^^, respectively. Also, the tree-size 
distribution (n^) behaves according to ~ s~^. Here, j^y and are the correlation length 
exponents in the directions parallel and perpendicular to the interface, respectively. Also, r is 
a critical exponent. However, due to the interplay between the discrete scale invariance of the 
underlying fractal substrates and the dynamics of the growing process, all these power laws are 
modulated by logarithmic periodic oscillations. The fundamental scaling ratios, characteristic of 
these oscillations, can be linked to the (spatial) fundamental scaling ratio of the underlying fractal 
by means of relationships involving critical exponents. We argue that the interplay between the 
spatial discrete scale invariance of the fractal substrate and the dynamics of the physical process 
occurring in those media is a quite general phenomenon that leads to the observation of logarithmic- 
periodic modulations of physical observables. 

PACS numbers: 68.35. Ct, 05.45.Df, 02.50.-r, 81.15.Aa 
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I. INTRODUCTION 



The study, characterization, and understanding of growth processes that take place un- 
der far-from equihbrium conditions are topics that have attracted great attention due to 
their relevance in many fields of science and technology s], 0, [4 1 . Growing aggregates of 
biological origin (bacteria, fungi, tumors, etc.), the deposition of thin films, the growth of 
magnetic materials, alloys and polymers, among others, have recently become the subject 
of extensive stud.es 31 Also, the dynamic evolution of .nterfaces is closely related 

to almost all growth processes [1]. The characterization of the properties of interfaces has 
achieved considerable progress during the last two decades, mostly due to the success of the 
concepts of the dynamic scaling theory developed by Family and Vicsek jo], lo|, which allows 
a comprehensive description of interfaces in terms of universality classes that group systems 
described by the same set of physically meaningful exponents. 

In order to provide a more complete description of growing aggregates, it is desirable not 
only to focus on the properties of the evolving interface, but also to attempt the simultaneous 
characterization of bulk properties. In fact, in some cases, growing processes lead to the 
formation of porous materials that inherently have very interesting physical and chemical 
properties with many potential practical applications. 

Within this broad context, ballistic deposition (BD), which was originally proposed by 



Void [llj as a model for the description of sedimentary rock formation, has become an 
archetypical system for the study of growing aggregates jl]. BD aggregates are character- 
ized by a porous structure in the bulk and a rough evolving interface. The interface rough- 
ness of BD has extensively been studied by means of computer simulations and analytical 
approaches in connection to the Kardar-Parisi-Zhang (KPZ) theory 

In spite of the considerable effort devoted to understanding growing aggregates in Eu- 
clidean substrates 

mm 

, to the best of our knowledge little attention has been drawn 
to the study of deposition models on fractal media. It should be expected that the interplay 
between the self-similarity of the substrates and the growing mechanisms would lead to the 
formation of interesting and complex porous (bulk) structures. Also, the fractality of the 
substrate would affect the self-affine nature of the growing interface. So, it would not be 
surprising to find nonvanishing scaling corrections to the well-established phenomenological 



dynamic approach of Family- Vicsek 
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10| . In fact, in some pioneering works [ij Q 
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and in more recently ones 
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171] it has shown that the power-law behavior of some ob- 



servables becomes modulated by logarithmic oscillations due to the effect of the underlying 
lattice structure. Also, in the last years two interesting books have been published where 



this issue is addressed 
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Within this context, the aim of this paper is to study the scaling behavior of the BD 
model on deterministic fractal substrates. The study is based on Monte Carlo numerical 
simulations analyzed by means of the Family- Vicsek phenomenological scaling approach 



0, Q 



9|, Il0| and on extensive numerical investigations of the scaling behavior of the internal 
structures of the growing system 2^, 21, 22, 23, 24, 25]. For this purpose the manuscript 



is organized as follows: firstly, in Section II we provide a brief theoretical background on 
the scaling behavior of growing aggregates. The fractal substrates used for the growth of 
BD aggregates and the concepts of space and time discrete scale invariance are described 
and discussed in Section III. Subsequently, in Section IV, we describe the BD model on 
fractal media, as well as the determination of the internal structure of the aggregates in 
terms of frozen and growing trees. The results obtained by applying the standard dynamic 
scaling approach to the data are presented and discussed in Section V, while the analysis of 
the data of the internal structure of the aggregates is performed in Section VI. Finally, our 
conclusions are stated in Section VII. 



II. BRIEF THEORETICAL BACKGROUND 

Models aimed to describe growing aggregates may be defined and studied by means of 
both continuous approaches, which involve the formulation of analytical equations, and 
discrete lattice models, which consider the deposition of individual particles. A discrete 
model is defined by means of a set of deposition rules that provides a detailed microscopic 
description of the evolution of the aggregate. In these discrete models, the growing interface 
of the aggregate is described by a discrete set h{i, t), which represents the height of site i at 
time t. The interface then has L'^ sites, where L is the linear size and d is the dimensionality 
of the substrate (as usual, d is assumed to be an integer). 

The dynamic evolution of the aggregate interface is characterized through the scaling 
behavior of the interface width W{L,t), given by 
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W{L,t) 



\ 



llL'Y}h(;i,t)-(h(tW. 



i=l 



For this purpose, the Family- Vicsek phenomenological scahng approach 
proved to be very successful, can be written as 



(1) 

which has 



iy(L,t) = L"iy*(t/L"), (2) 

where W* is a scaling function. In fact, it may be expected that W{L,t) will show the 
spatiotemporal scaling behavior given by [91, ll0|: W (x for t > tc and W{t) oc for 
t -C tci where t^. oc is the crossover time between these two regimes. The scaling exponents 
a, /3, and z = a/ (3 are called roughness, growth and dynamic exponents, respectively. Also, 
different models can be grouped into universality classes when they share the same scaling 
exponents. 

An alternative method that can also be used for the characterization of interfaces is 
related to the description of the internal structure of the growing system. This approach is 
based on the fact that any growing system can effectively be rationalized on the basis of a 
treeing process, i.e., a ny g rowing structure can be thought as the superposition of individual 



trees 120 
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26| . Those trees that spread out incorporating additional growing 



centers, e.g., capturing particles, developing new branches, are said to be alive. In contrast, 
other trees that may stop growing due to shadowing by surrounding growing trees are termed 
dead trees. The structure of dead trees remains frozen because it cannot be modified by any 
further growth. It is well known that for growing aggregates on substrates having integer 
dimension, the rms height ( hg ) and the rms width { Wg) oi dead trees of size s [s is the 
number of particles belonging to the tree) obey simple power laws given by 



K oc s'^ii 



(3) 



and 



Wg (X s 



(4) 



where v\\ and z/^ are the correlation 
growing direction of the aggregate 



length exponents parallel and perpendicular to the main 
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26|, respectively. 



By assuming both the proportionahty between the correlation length perpendicular to 
the main growing direction and Ws, as well as hs oc t, one has that the dynamic exponent 
z = a/l3 is given by 2^] 



Z/||/Z/_L. 



(5) 



Furthermore, one also expects that during the competition among trees along the evo- 
lution of the aggregate, the existence of large neighboring trees may inhibit the growing 
of smaller ones. This competing process ultimately leads to the death of some trees that 
become frozen within the underlying aggregate. These prevailing large trees continue the 
competition within more distant trees in a dynamic process. Since this situation takes place 
on all scales, it is reasonable to expect that the tree size distribution (n^) should also exhibit 



a power-law behavior, so that 
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ris^ s ^, (6) 

where r is an exponent. 

III. THE FRACTAL SUBSTRATES AND DISCRETE SCALE INVARIANCE 

In the present paper we used Sierpinski carpets (SC) as substrates for the growth of bal- 
listic aggregates. In fact, SC's provide generic models for the building of both deterministic 
and nondeterministic fractals. In order to generate a SC embedded in c? = 2 dimensions, one 
proceeds as follows: a square is divided into subsquares, and then {P — Nocc) subsquares 
are deleted from the initial square {Nocc is the number of occupied subsquares). This pro- 
cess is iterated in the remaining subsquares k times, where k accounts for the number of 
different generations. If the deleted subsquares are chosen in the same way in all iterations, 
the resulting fractal is deterministic, but if the deleted subsquares are selected at random, 
one generates a nondeterministic fractal. The mathematical fractal, obtained in the limit 
k oo, is generically called SCx{l, Nocc)- Also, the fractal associated with a finite number 
of segmentation steps is denoted by SC^i^l, Nocc, k). In both cases the index x refers to the 
topological features of the generating cell {SCx{l, Nocc, !))• The size of the lattice, where the 
finite fractal is embedded, then is L = /'^. 

5 



For deposition models grown in integer-dimensional substrates one always has that all 
fragments or parts of the aggregate are connected to each other through paths of nearest- 
neighbor occupied sites. So, it is no longer possible to have isolated fragments of the ag- 
gregate on the substrate. On the other hand, for some fractal substrates with noninteger 
fractal dimension, it could be possible to observe the formation of isolated fragments of the 
aggregate, and consequently in order to avoid this shortcoming, one has to carefully select 
suitable substrates. So, this point is essential for the choice of the fractals that can be used 
in order to study ballistic deposition on this kind of substrate. Accordingly, in this paper 
we use deterministic fractal substrates generated by taking / = 3 and Nocc = 5, and in order 
to prevent the fragmentation of the aggregates, as mentioned above, and to account for the 
usual requirement of periodic boundary conditions, one ends up with only three generating 
cells, as shown in figure [H 




FIG. 1: Sketches of the three different generating cells used to grow ballistic aggregates by keeping 
/ = 3 and Nqcc = 5. 

Let us now discuss on the concepts of fractal dimensionality and discrete scale invariance 
associated with the SC's used. The fractal dimension df characterizes the dependence of 
the mass M{L) (or equivalently, the number of occupied sites of the fractal) as a function 
of the linear size L of the system, so that if we now consider an amplification of the system 
of size bL, one has 
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M{bL) = b'^fM{L), 

which yields the following solution 



(7) 



M{L)=BL'^f, (8) 

where i? is a constant and df is a noninteger dimension known, after Mandelbrot, as the 
fractal dimension. By applying equation (JTj) to the carpets SCx{l, Nocc, k) one gets df = 
ln{Nocc) / in{l) , which leads to df = ln{5)/ln{3) ^ 1.465 for the carpets shown in figured] 
with I = 3 and Nocc = 5. 

In general, the factor b in equation ([7]) could be an arbitrary real number, leading to 
continuous scale invariance. However, deterministic fractals exhibit discrete scale invariance 



(DSI) 28|, which is a weak kind of scale invariance such that b is no longer an arbitrary real 
number, but it can only take specific discrete values of the form 6„ = where bi is a 

fundamental scaling ratio. Then, for the case of DSI, the solution of equation ([7]) yields 

M{L) = L'fF(p^], (9) 



Jog(^i), 

where F is a periodic function of period one. Notice that for the SCx{l, Nocc, k) one has that 



bi = I. The measurement of soft oscillations in spatial domain 28| is a signature of spatial 
DSI. 

Very recently one of us found evidence of discrete scale invariance in the time domain 
)y measuring the relaxation of the magnetization in the Ising model on Sierpinski carpets 
Subsequently, it has been conjectured that physical processes characterized by 
an observable 0{t), occurring in fractal media with DSI, and that develop a monotonically 



increasing time-dependent characteristic length ^(t), may also exhibit time DSI 17|]. In fact, 
by assuming 

eoct^/"°, (10) 
where zd is a dynamic exponent, it can be shown that 0{t) has to obey time DSI according 



to 
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O(t)^Ct^/-F(0 + Jg^), (11) 
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where C and are constants, and 7 is the relevant exponent in the expected power-law 
behavior of the observable 0{t). So, the conjecture given by equation (flTl) implies the 
existence of a logarithmic periodic modulation of time observables characterized by a time- 
scaling ratio T given by 



T = 6f, (12) 

see also equation (Q. It is worth mentioning that in the case of growth models, equation 
ffTOj) can be identified, e.g., with the time development of the correlation length along the 
direction parallel to the interface given by oc t^^^. 



As follows from figure [T], in this work only fractals with finite ramification order 30|] are 
used. A finite ramification implies that the fractal structure has weak points, where only a 
finite number of links connect two parts of arbitrary size together. 



IV. THE BD MODEL ON FRACTAL SUBSTRATES AND DEFINITION OF ITS 
INTERNAL STRUCTURE 

The lattice version of BD is simple to describe: particles fall vertically onto the substrate 
from a random position above the surface. When a particle reaches the surface, it sticks on 
the first site encountered that is a nearest-neighbor of an already deposited particle. Due to 
this constraint the growth of an interface essentially parallel to the substrate is observed. For 
substrates of integer dimension, the BD model can be described by the Kardar-Parisi- Zhang 



(KPZ) equation l|, ll2|, namely, 



= ^ + z/,V^/i(x, t) + ^[VMx, t)]' + r7(x, t). (13) 

In this equation the nonlinear term represents the lateral growth or the appearance of a 
driven force, Uo accounts for the effective surface tension, JF is the flux of incoming particles, 
and ?7(x, t) is a Gaussian noise with zero configurational average. In d = 1 dimension, 
equation f|T3l) can be solved exactly and the resulting exponents are z = 3/2, a = 1/2 and 
f3 = 1/3. For d = 2 dimensions, one still lacks an analytical solution, but according to 
numerical simulations one has estimations of the exponents given by a ~ 0.40 and j3 ~ 0.24 
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In order to describe the internal structure of BD aggregates, one considers that trees are 
formed by assuming that any newly deposited particle belongs to the same tree as that of the 
nearest neighbor particle where it is attached 3]. Also, if the deposited particle has more 
than one nearest neighbor belonging to different trees, one of them is selected at random 
and the particle is incorporated into that tree. Relevant exponents (see equations ([3]), (jl]), 
and ([6])) have already been determined yielding i/y = 0.60(1), z/_l = 0.40(1), and r = 1.40(1) 
in d = 1 dimension, and uh = 0.45(1), z/_l = 0.29(1), and r = 1.57(1) in d = 2 dimensions 



[see, e.g. 



24| and references therein). 



On the other hand, for BD growth on fractal substrates, particles are deposited randomly 
on occupied sites of the fractal only by following the rules described above for the case 
of integer-dimensional substrates. Also the building procedure employed to determine the 
trees used for the description of the internal structure of the whole aggregate is independent 
of the substrate. 

For the purpose of a numerical simulation, the Monte Carlo time step (mcs) involves the 
deposition of L'^f particles. 



V. APPLICATION OF THE FAMILY- VICSEK PHENOMENOLOGICAL DY- 
NAMIC SCALING APPROACH 

Figure [2] shows log- log plots of W versus t obtained for the BD model on fractal substrates 
obtained by using the generating cell SCa{3, 5, 1). For short times, say t < 3 mcs, the random 
growth of the interface is observed because the random deposition (RD) process dominates. 
At this stage, correlations have not been developed yet and one has that, according to 
equation ([2]), W(t) oc t^^-", with Prd = 1/2, holds. During an intermediate time regime, say 
3 mcs < t < tc, correlations develop since the BD process now prevails, leading to the typical 
growth regime W{t) oc t^^°, see also equation ([2]). At a later stage, for t > tc, correlations 
can no longer develop due to the geometrical constraint of the lattice size, and the saturation 
regime {Wsat{L) oc L^-^o) is observed, as expected from equation ([2]). Here, tc oc is the 
crossover time between the growing and the saturation regimes of the interface width. 

In order to obtain the roughness, growth, and dynamic exponents, we followed the stan- 
dard procedure [Ij], e.g., the inset in figure [2] shows a log-log plot of WsatiL) versus L. Here, 
a slight but noticeable systematic upward deviation of the data is observed. However, an 
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FIG. 2: Log-log plots of the interface width (W) versus time for the BD model on SCai^"^ ,5,k) 
with k = 2 — 6. The inset shows a log-log plot of Wsat{L) versus L. The dotted line with slope 
0.47 corresponds to the best fit of the data. 

effective exponent (a) can be determined and the best fit of the data yields a = 0.47 ±0.04. 
Also, we determined (3 = 0.28 ± 0.03 and z = 1.56 ± 0.05 for the BD model on the fractal 
substrate obtained by using the SCa{S, 5, 1) generating cell. 

Figures [3] and H] show log-log plots of W versus t obtained for the BD model on fractal 
substrates built up by using the generating cells SCh{3,5, 1) and S'Cc(3,5, 1), respectively. 
Again, the standard procedure was attempted in order to determine the roughness, growth, 
and dynamic exponents. Accordingly, the insets in figures [3] and H] show log-log plots of 
Wsat(yL) versus L. For the case of the SCb{S,5,l), again, a noticeable systematic upward 
deviation of the data is observed, but an effective exponent (a) can be determined, yielding 
a = 0.44 ± 0.04. Also, we determined f3 = 0.26 ± 0.03 and z = 1.55 ± 0.05. On the 
other hand, for the SCc{3, 5, 1) we observed strong systematic deviations of the data that 
prevent the evaluation of the exponents. This results strongly suggest that in the case of 
aggregates grown on fractals generated by SCc{3, 5, 1) either the finite-size effect are very 
important and one has to introduce important scaling corrections to the standard Family- 
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FIG. 3: Log-log plots of the interface width (W) versus time for the BD model on SCi,{3'^,5,k) 
with k = 2 — Q. The inset shows a log-log plot of WgatiL) versus L. The dashed line with slope 
0.44 corresponds to the best fit of the data. 

Vicsck approach, or (instead of we are not be able to identify two competing processes leading 
to the formation of the aggregate) eventually the onset of a crossover effect is present. The 
understanding/clarification of which one of the two possibilities may be true is a very difficult 
task, which deserves further work, and is beyond the scope of the present paper. 

Anyway, the obtained effective exponents are almost the same (within error bars) for both 
the SCai3, 5, 1) and the SCbiS, 5, 1), namely a a; 0.46, /? w 0.27, and z « 1.55. It is worth 
mentioning that these exponents nicely interpolates between the exact values corresponding 
to d = 1 and the best estimates reported for d = 2, namely 0.24((i = 2) < /3 ~ 0.27 < 
l/3{d = 1), 0.40(ci = 2) < a 0.46 < l/2{d = 1), and 1.67(ci ^ 2) > z ^ 1.55 > 3/2{d = 1). 
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FIG. 4: Log-log plots of the interface width (W) versus time for the BD model on S'Cc(3^,5,fc) 
with k = 2 — 6. The inset shows a log-log plot of WsatiL) versus L. The dashed line with slope 
0.48 have been drawn for the sake of comparison. 

VI. CHARACTERIZATION OF THE INTERNAL STRUCTURE OF THE AG- 
GREGATES 

The dependence of the rms width and rms height of the trees forming the aggregates 
on the tree size s, for the cases of the SCb{S,5, k) substrates {k = 3,4,5), is shown in log- 
log plots in figures [5] and [6l respectively. For the growth of the BD model on^nonfractal 
substrates it is known that the power laws given by equations (I3l) and (jlj) hold [24], but in 
the case of fractal media, one also clearly observes soft oscillations with a logarithmic period, 
which modulate the power laws. Figures [5] and [6] also show a direct relationship between 
the number of observed oscillations and the generation {k) of the Sierpinski carpets used as 
substrates. 

Figure [7] shows log-log plots of the tree-size distribution functions corresponding to the 
BD model grown on SCb{3,5, k) substrates {k = 2,3,4,5). Again, for the BD model on 
nonfractal substrates the power law given by equation ([6]) holds, but when a fractal carpet 
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FIG. 5: Log- log plots of rms width of the trees as a function of the tree size s for the BD model 
on 5Cfe(3, 5,/c) substrates {k = 3,4,5). 

is used as substrate, a soft oscillation with a logarithmic period modulates the power-law 
behavior. Also, there is a one to one relationship between the number of observed oscillations 
and the number of generations of the carpet. 

After a systematic study, we found that the rms width, the rms height, and the tree-size 
distribution for BD aggregates grown on different deterministic fractal substrates, generated 
by / = 3 and Nocc = 5, exhibit a quite similar behavior, namely, soft oscillations with a 
logarithmic period, which modulate the expected power laws. Figures [HI El and fTOl show log- 
log plots of the rms width, the rms height, and the tree-size distribution, obtained for these 
fractal substrates, respectively. Table I summarizes the exponents obtained by fitting the 
modulated power laws describing the properties of the internal structure of the aggregates. 
Also, data corresponding to the exponent z obtained by means of the Family- Vicsek dynamic 
scaling approach have been included for the sake of comparison. The values of the exponents 
F = z/||, i'±, and r, have been obtained by fitting the data with the function 
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FIG. 6: Log- log plots of rms height of the trees as a function of the tree size s for the BD model 
on 5Cb(3, 5, k) substrates {k = 3, 4, 5). 



= Y + J2bnCos{nu;Hs) + e„), (14) 

^ ' n=lr 



where X = hg, Wg, and n^, respectively [31j- Here, fe„ and ^„ are constants. When the 
fit is performed by taking > 1, the presence of higher harmonics to the fundamental 
frequency is considered. So, the values of the exponents obtained by using = 1 could 



be improved 



311 ]. However, when we considered higher harmonics, inspite of the fact that 



typically we have ~ 0.2, the values of the exponents obtained by fitting the data with 
= 1, and 2 are almost indistinguishable (within errors bars). 



In order to rationalize our findings, let us discuss the soft oscillations observed in the rms 
width mentioned above. Here, we can identify two logarithmic periods, one for the size of 
the trees (As) and the other for the rms width (A^J, as is schematically shown in figure [TTl 
These periods become evident by calculating the derivative of ln{ws) with respect to ln(s), 
and the derivative of ln(s) with respect to ln{ws), as shown in figure [12] for the case of the 
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FIG. 7: Log-log plots of the tree size distribution corresponding to the BD model on 5C;,(3, 5,A:) 
substrates (A; = 2, 3, 4, 5) 



TABLE 1: List of exponents describing the properties of the internal structure of ballistic aggregates 
on different fractals. The exponent z (5*^ column) is obtained by using equation ^ and the values 
of vx_ and are listed in the 2"^*^ and 3'''^ columns, respectively. Also, the exponent z column) 
is evaluated by using data obtained by applying the Family- Vicsek dynamic scaling (equation ([2])). 
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z = a/ (3 


5Ca(3,5,l) 


0.34(1) 


0.54(1) 


1.48(1) 


1.59(6) 


1.56(5) 


5^(3,5,1) 


0.34(1) 


0.55(1) 


1.48(1) 


1.60(6) 


1.55(5) 


5Ce(3,5,l) 


0.35(1) 


0.56(1) 


1.48(5) 


1.61(6) 





SC},{?>, 5, 5) substrate. Also, it is worth mentioning that by following the same procedure, 
we can identify one logarithmic period for the rms height (A/^J and another one for the 
size distribution (A„J. Furthermore, as mentioned above, there is a one to one relationship 
between the number of observed oscillations and the generations of the carpets. So, as all 
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FIG. 8: Log- log plots of rms width of the trees as a function of the tree size s for the BD model 
on 5Cx(3,5,5) substrates {x = a,b,c). 

measurements share the same range of values of s, the period of the oscillation observed 
for the size of the trees (A^), when k tends towards infinity, has to be the same for all 
observables, namely, the rms width, the rms height, and the tree-size distribution. 

Now, based on both the already discussed results and equation (jl]), we conjecture that 
for the case of a fractal substrate and in the limit k oo, the derivative of ln(w5) with 
respect to ln(s) is a periodic function (with a logarithmic period) that modulates a constant 
value. So, the conjecture can be written as 



dhi{ws 
d\n{s) 



z/x + C(ln(s)), 



(15) 



where C,{x) is a periodic function. By integration of equation f[T5]) it follows that 



= As^^exp'^*(i^(^)), (16) 

where A is a constant and Ci.^)* is a periodic function such that dC,* {\n.{s)) / dhi^s) = ^(ln(s)). 

Let us now take two tree sizes, Sa and Sb, separated by a period (on a logarithmic 
scale), as shown in the horizontal axis of figure [TTl Also, their corresponding rms widths, 
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FIG. 9: Log- log plots of rms height of the trees as a function of the tree size s for the BD model 
on SCx{3, 5, 5) substrates {x = a, b, c). 

called Wa and Wb, respectively, are separated on a logarithmic scale by a period A^,^ (see 
figure fTTl) ■ Then, by assuming that equation (|T6l) holds, it can be found that 

A,„, = i^±Xs- (17) 

The existence of a logarithmic period in the rms width of the trees (A^J can be understood 
due to the fact that each tree can only spread over the fractal and, consequently, the rms 
width of the trees is constrained by the geometrical features of the underlying structure. In 
the case of a deterministic Siernspisky carpet, that structure is constructed by the iteration 
of a generating cell, so that the topological details of the generating cell are present in all 
the sample. In our simulations, a generating cell of side / = 3 leads to the occurrence of 
discrete scale invariance, such that the underlying structures are self-similar only at scales 
of size Z", where n is an integer. Based on theses concepts, it is expected that A^^ = In (3). 
In order to check this statement, in the inset shown on the left-hand side of figure [I2] we 
have drawn a line of size ln(3), which nicely fits the logarithmic period of the oscillation 
of the rms width. This result has also been carefully checked for all the studied cases by 
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FIG. 10: Log- log plots of the tree size distribution corresponding to the BD model on SCx{3, 5, 5) 
substrates {x = a, b, c). 



measuring the peak-to-peak distance of the modulating oscillations. This finding leads us 
to conjecture that, for aggregate grown on deterministic fractal substrates generated by a 
generating cell SC^il, Nocc, 1), the logarithmic period in the rms width of the trees should 
obey the following relationship 



A^. =ln(/). (18) 

Equation (fT8l) refiects the fact that the size of the trees in the direction parallel to the 
substrate exhibits spatial discrete scale invariance with the same fundamental scaling ratio 
as that of the underlying fractal, namely, bi = I (see also equation ([9j) ). 

In order to check this conjecture, figure [13] shows a log-log plot of the rms width as a 
function of the tree size s for the case of an SCb{2, 3, 8) substrate. In this case, the size of 
the generating cell is / = 2, so according to the conjecture we expect A^;^ = ln(2). The inset 
placed on the left of figure [13] shows the derivative of ln(s) with respect to \ia{ws), and a fine 
of size ln(2), which nicely fits the value expected for A^^, has been drawn. 

In order to understand the oscillatory modulation of the power laws observed in the rms 
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FIG. 11: Schematic view of the dependence of the rms width of the trees {wg) versus the tree size 
(s), showing both the logarithmic period of the size of the trees (A^), and the logarithmic period 
of the rms width (A^^ ) . 

height of the trees, we followed the same method already discussed for the case of the rms 
width. Figure [H] shows the rms height as a function of the tree size s, as obtained for 
the SCb{3, 5, 5) substrate. Also, both the derivative of lia{hs) with respect to ln(s) and 
the derivative of ln(s) with respect to ln{hs) are shown in the insets placed on the lower 
right-hand and the upper left-hand sides of figure [HI respectively 

So, based on equation ([3]) and the obtained results, we now conjecture that 



d\n{hs 
d\n{s) 



v\\ +ri{\n{s)), 



(19) 



where ri{x) is a periodic function. Then, after integration one gets 



/i, = Ss"iiexp^*('"(^)), (20) 

where S is a constant and ri{x)* is a periodic function such that 5?7*(ln(s))/91n(s) = ?7(ln(s)). 

By using equation (!20l) and following the method already applied to the case of the rms 
width, it can also be found that 
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FIG. 12: Log- log plots of the rms width of the trees as a function of the tree size s as obtained 
for balhstic aggregates on 5Cb(3^,5,5) substrates. The derivative of ln(u;s) with respect to ln(s) 
is shown in the inset placed at the lower-right hand side, while the derivative of ln(s) with respect 
to ln{ws) is shown in the inset placed on the upper-left hand side. The dashed line in the inset 
placed at the lower-right hand side correspond to the best fit of the data using Eq. [TH 



Xh. = iy\\Xs. (21) 
Then, by using equations f|T7|) . ffTSl) . and fl2T]) one obtains 

A;,, = ^ ln(3) = z ln(3) = z ln(/) = z (22) 

where equation ([5]) has also been used. Let us now recall that, according to our conjecture 
given by equation ( fTTl) . it should be expected that the rms high of the trees, which is an ob- 
servable that evolves along the growing (time) direction, will be coupled to the space discrete 
scale invariance of the underlying fractal, then exhibiting time discrete scale invariance. In 
fact, the coupling is nicely reflected by equation (!22|) . which is the realization of equation 
( [T2|) . with T = exp(A/iJ being the fundamental time scaling ratio. 



20 



10° 10^ 10^ 10^ 10' 

s 



FIG. 13: Log-log plots of the rms width as a function of the tree size s as obtained for the 
S'Cf,(2,3, 9) substrate. The derivative of ln(s) with respect to ln{ws) is shown in the inset placed 
on the left. 

Finally, the same procedure can be applied to understand the oscillation observed around 
the power laws of the size distribution. In fact, figure [I5] shows the tree size distribution 
function for the SCb{3, 5, 5) substrate. The derivative of In(ns) with respect to ln(s) is shown 
in the inset placed on the left-hand side, while the derivative of ln(s) with respect to In(ns) 
is shown in the inset placed on the right-hand side. 

Based on the obtained results, now we conjecture that in the case of a fractal substrate 
where k tends towards infinity, the logarithmic derivative of the size distribution with respect 
to ln(s) should be a periodic function around a constant value. So, in this case one has 

^ (-) 

where 6{x) is a periodic function. Then, after integration one gets 

n, = Cs"exp'*(i"(^», (24) 
where C is a constant and 9{x)* is a periodic function such that d9* {ln{s)) / dlia{s) = 6'(ln(s)). 
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FIG. 14: Log- log plots of the rms height as a function of the tree size s as obtained for the 
SCh{3,5,5) substrate. The derivative of ln{hs) with respect to ln{s) is shown in the inset placed 
on the right, and the derivative of the ln(s) with respect to ln(/is) is shown in the inset placed on 
the left. 

Now, by using equation (!24l) and following the same method as in the case of the rms width, 
it is easy to find that 

A„, = tXs. (25) 

Table II summarizes the values of the different logarithmic periods obtained for the BD 
model on SCx{'i, 5, 5) substrates (x = a, 6, c). The results were obtained by fitting the data 
shown in the insets of figures 12, 14 and 15, with the aid of equation [14] for N = 2. In 
the case of the BD model on SCa{3, 5, 5) fractals, the data exhibit some fluctuations and it 
is no longer possible to determine a reliable value of A^^. The same shortcoming is found 
when attempting the evaluation of Xh^ for the case of the SCc{3, 5, 5) substrate. In order 
to overcome this shortcoming, it should be necessary to perform the simulations on fractals 
with k ^ 6, but this task is very CPU-time demanding, so it is beyond our computational 
capabilities. 
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FIG. 15: Log-log plot of the tree size distribution as obtained for the 5(7^(3,5,5) substrate. The 
derivative of In(ns) with respect to In(s) is shown in the inset placed on the right, and the derivative 
of ln(s) with respect to In(ns) is shown in the inset placed on the left. 



TABLE 11: List of the logarithmic periods obtained by fitting the log-periodic modulations of the 
power laws describing the internal structure of ballistic aggregates on different fractals. The data 
are also compared to scaling relationships that are explained in the text. 



Fractal 




ln(/) 


Xs 






z\n{l) 




Xns 


tXs 


SCa{3,5,5) 




1.0986 


3.4(1) 




1.80(8) 


1.75(6) 


1.84(6) 


4.8(2) 


5.0(2) 


SCt{3,5,5) 


1.10(3) 


1.0986 


3.2(1) 


3.2(1) 


1.70(8) 


1.76(6) 


1.77(6) 


4.7(2) 


4.7(2) 


5Cc(3,5,5) 


1.10(3) 


1.0986 


3.3(1) 


3.2(1) 




1.77(6) 


1.84(6) 


4.7(2) 


4.9(2) 



It is worth mentioning that all determined logarithmic periods can be compared to scaling 
relationships also involving exponents determined independently. In fact, the values of A^^ 
(2"*^ column) are compared to ln(/) (3''°' column), as it follows from equation f|T8|) . obtaining 
an excellent agreement. Also, according to equation (fTTll . the period (4*^^ column) can 
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be compared with the quotient given by Xws/^i- (S*'* column), where the exponent z/_l is 
taken from Table I, leading to an excellent agreement. On the other hand, the period 
(6*^ column) can be compared with the relationships given by equations (1221) (7*'' column) 
and fl^Tl) (S*'' column), where the exponents z and v\\ are taken from Table I, obtaining 
agreement within error bars. Finally, the relationship X^^^ = tX^ given by equation (1251) . 
can be verified by comparing columns 9^^ and 10*'', also observing agreement within error 
bars. 

Finally, we would like to comment that for growing aggregates in (i— dimensional homo- 
geneous media, one expects that the average volume of the frozen trees of size s (vs), will 
scale according to [24 1 



V, oc oc s'^ii+^'^-i)'^^, (26) 

so that one can define f ^ oc s'^, where vr is an exponent. In previous work we showed that for 
ballistic aggregates in 1 < d < 5 one has that vr ~ 1, indicating that the trees are compact 
structures [2J]. By using this procedure, we determined vr = 1.05(2), 1.06(2), and 1.07(2) for 
the fractals of type a, b, and c (see figured]), respectively. Since we evaluate the statistical 
error only, these exponents could be considered as preliminary evidence that frozen trees are 
also compact for ballistic aggregates in fractal media. 



VII. CONCLUSIONS 



We studied the ballistic deposition growth model in deterministic fractal substrates, em- 
bedded in (i = 2 dimensions, by means of numerical simulations, analyzing the data according 
to well established scaling theories. 

While the standard Family- Vicsek scaling approach seems to hold at least qualitatively, 
i.e., the interface width exhibits an initial power-law increase followed by saturation, sys- 
tematic deviations from the expectation Wsat{L) oc L" are observed for large L— values. In 
fact, the measured values of the interface width are larger than the scaling predictions, sug- 
gesting that the standard scaling approach has to be modified in order to properly describe 
the behavior of the interface upon growth on fractal media. Anyway, the obtained (effective) 
scaling exponents a, (3 and z, interpolate between the already known values corresponding 
to d = 1 and d = 2 Euclidean dimensions. Furthermore, the relationship z = a/(3 holds, 
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within error bars, and the determined dynamic exponent z is in reasonable agreement with 
subsequent independent measurements, performed by studying the internal structure of the 
aggregates (see, e.g.. Table I). 

On the other hand, the study of the properties of the internal structure of the aggregates, 
in terms of scaling laws describing the properties of frozen trees, i.e., branched structures that 
are inhibited for further growth due to shadowing effects and lie below the interface, provides 
an interesting and rich physical behavior. In fact, the standard scaling theory applied 
to the deposition on substrates of integer dimension predicts the power-law behavior of 
relevant physical observables as a function of the size of the trees, such as the rms width and 
height of the trees, as well as for the tree size distribution. However, for fractal media with 
dj = ln(5)/ln(3) and dj = ln(3)/ln(2), we show that these power laws become modulated 
by soft log-periodic oscillations. These oscillations reflect the interplay between the spatial 
discrete scale invariance (DSI) of the underlying fractal substrates and the growing process. 
While it is well known that dynamic observables become coupled to the spatial DSI given 



time DSI 
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171], growing aggregates in fractal media provide a richer scenario. On the 



one hand, the rms width of the frozen trees, i.e., an intrinsically spatial observable measured 
along the direction parallel to the substrate, exhibits spatial DSI with the same fundamental 
scaling ratio (In A^^ = ln(/) ) as that of the fractal. On the other hand, the rms height 
of the frozen trees, which is measured along the growing (time) direction perpendicular to 
the substrate, exhibits time DSI and its fundamental scaling ratio (In X^^ ) is coupled to 
InAu;^ through the dynamic exponent z (see equation ( l22i) ). which is the realization of the 
more general relationship, also valid for any dynamic observable, given by equation (|TOl) . 
Of course, the size distribution of the frozen trees also exhibits DSI, which in a subtle way 
involves both spatial and time DSI, since it is an observable that depends on the structure 
of the trees developed along the directions parallel and perpendicular to the substrate. 

We hope that this work will contribute to the understanding of growing aggregates in 
fractal media, also stimulating theoretical work aimed to describe the interplay between the 
underlying symmetries of the substrates and the behavior of physical observables character- 
istic of the processes occurring in those media. 
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